Cancellation exponent and multifractal structure in two-dimensional magnetohydrodynamics: 
direct numerical simulations and Lagrangian averaged modeling 
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We present direct numerical simulations and Lagrangian averaged (also known as a-model) simulations of 
forced and free decaying magnetohydrodynamic turbulence in two dimensions. The statistics of sign cancel- 
lations of the current at small scales is studied using both the cancellation exponent and the fractal dimension 
of the structures. The a-model is found to have the same scaling behavior between positive and negative con- 
tributions as the direct numerical simulations. The a model is also able to reproduce the time evolution of 
these quantities in free decaying turbulence. At large Reynolds numbers, an independence of the cancellation 
exponent with the Reynolds numbers is observed. 
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The magnetohydrodynamic (MHD) approximation is of- 
ten used to model plasmas or conducting fluids in astrophys- 
ical and geophysical environments. However, given the huge 
amount of temporal and spatial scales involved in the dynam- 
ics of these objects, simulations are always carried out in a 
region of parameter space far from the observed values. La- 
grangian averaged magnetohydrodynamics (LAMHD), also 
called the MHD alpha-model Ql |21 (or the Camassa-Holm 
equations in early papers studying the hydrodynamic case 
131), has been recently introduced as a way to reduce the 
number of degrees of freedom of the system, while keeping 
accurate evolution for the large scales. This approach (as 
well as large eddy simulations, or LES, for MHD; see e.g. 
01) is intended to model astrophysical or geophysical flows 
at high Reynolds numbers using available computational re- 
sources. Several aspects of the MHD alpha-model have al- 
ready been tested in two and three dimensions at moderate 
Reynolds numbers, against direct numerical simulations of the 
MHD equations |2]. These studies were focused on compar- 
isons of the evolution of global quantities and the dynamics of 
the large scale components of the energy spectrum Q,|5|l. 

All these models introduce changes in the small scales in 
order to preserve the evolution of the large scales. In several 
cases, it is of interest to know the statistics of the small scales. 
It is also important to model properly the small scales because 
they have an effect on large scales, as for example in the case 
of eddy noise: the beating of two small scales eddies produces 
energy at the large scale, and this may affect the global long- 
time evolution of the flow, an issue that arises in global climate 
evolution or in solar-terrestrial interactions. Moreover, plas- 
mas and conducting fluids generate thin and intense current 
sheets where magnetic reconnection takes place. In these re- 
gions, the magnetic field and the current rapidly change sign, 
and after reconnection the magnetic energy is turned into me- 
chanical and thermal energy. These events are known to take 
place in the magnetopause 01, the magnetotail fl^, the solar 
atmosphere |^, and the interplanetary medium ||^. 

Current sheets are strongly localized and intermittent. To 
preserve reliable statistics of these events in models of MHD 
turbulence is of utmost importance to model some of these as- 
trophysical and geophysical problems. In this work, we study 
whether the MHD alpha-model is able to reproduce the statis- 



tics and scaling observed in these phenomena. 

In order to measure fast oscillations in sign of a field on ar- 
bitrary small scales, the cancellation exponent was introduced 
fTo , 11, 121. The exponent is a measure of sign-singularity. 
We can define the signed measure for the current (x) on a 
set Q{L) of size L as 



Q,(i) 



c^x|j^(x)| 



(1) 



QiL) 



where {Qi{l)} C Q{L) is a hierarchy of disjoint subsets of 
size I covering Q{L). The partition function x measures the 
cancellations at a given lengthscale I, 



x{i) 



Q^(l) 
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Note that for noninteger L/l the subsets will not cover Q{L) 
and finite size box effects must be considered in the normal- 
ization of Eq. Q. We can study the scaling behaviors of the 
cancellations defining the cancellation exponent k, where 



(3) 



Positive K indicates fast changes in sign on small scales (in 
practice, a cut-off is always present at the dissipation scale). 
A totally smooth field has k = {). This exponent can also be 
related with the fractal dimension D of the structures fl^ . 



K^{d~D)/2, 



(4) 



where d is the number of spatial dimensions of the system. In 
some circumstances, we will also be interested on the cancel- 
lation exponent for the vorticity cj^. In that case the vorticity 
replaces the current in the definition of [Eq. Q]. 

Under special assumptions, relations between the cancella- 
tion exponent and scaling exponents have also been derived 
fllll . Positive cancellation exponent k has been found in 
plasma experiments fl^ . direct simulations of MHD turbu- 
lence fl2il . in situ solar wind observations fisll . and solar pho- 
tospheric active regions fl4ll . where changes in the scaling 
were identified as preludes to flares. 

In this work we will consider both free decaying and forced 
simulations of incompressible MHD and LAMHD turbulence 



2 



in two dimensions (2D). The MHD equations in 2D can be 
written in terms of the stream function and the z component 
of the vector potential A^, 



dtA, 



i^V% (5) 
(6) 



where the velocity and magnetic field are given by v = 
V X ) and B = V x (^^^z) respectively, and \F, G] 
dxFdyG — dxGdyF is the standard Poisson bracket. The 
LAMHD equations are obtained by introducing a smooth- 
ing length a, and the relation between smoothed (denoted 
by a subindex s) and unsmoothed fields is given by F = 
(1 - v2)Fs, for any field F. The system of LAMHD equa- 
tions in this geometry j^] is 



dtAs^ 



j/V^'l' (7) 
(8) 



For both systems of equations, the current is given by jz = 
— V^^z, and the vorticity by lOz = — V^^*. In these equations 
and in all the following figures, all quantities will be given 
in familiar Alfvenic dimensionless units. Equations ( I5l8t are 
solved in a periodic box using a pseudospectral code as de- 
scribed in The code implements the 2/3-rule for dealias- 
ing, and the maximum wavenumber resolved is kmax = N / 3, 
where N is the linear resolution used in the simulation. All 
the fields are written in dimensionless units. 

To characterize the oscillating behavior and sign singular- 
ities in the flows obtained from the MHD and LAMHD sim- 
ulations, we perform a signed measure analysis and compute 
the cancellation exponent k for the current and for the vor- 
ticity. Following Eq. Q, its value is obtained by fitting 
x(l) = c(llV)~'^ through the inertial range, where L = 2Tr 
is the length of the box, and c is a constant. The lengthscales 
in the inertial range used for this fit are obtained studying the 
scaling of the third order structure function fl^ . 

We first present results for a forced MHD simulation with 



1024 grid points, with 77 



1.6 X 10" 



Both the mo- 



mentum and the vector potential equations were forced. The 
external forces had random phases in the Fourier ring between 
k = 1 and k = 2, and a correlation time of A< = 5 x 10^^. 
The system was evolved in time until reaching a turbulent 
steady state. The amplitude of the magnetic force averaged 
over space was held constant to 0.2, and the amplitude of the 
mechanical force to 0.45, in order to have the system close to 
equip artition. Two more simulations using the LAMHD sys- 
tem were carried out, with the same parameters as the MHD 
run but with resolutions of 512^ grid points {a w 0.0117), 
and 256^ grid points (a w 0.0234) respectively (the choice 
a — 2 /kmax is conventional Q,!!!). The Kolmogorov's ki- 
netic and magnetic dissipation wavenumbers in the MHD run 



are k^ 



332; in all the LAMHD simulations these 



wavenumbers are larger than the largest resolved wavenum- 
ber kmax, by virtue of the model. Note that although it is 
common to reduce the spatial resolution even more in stud- 
ies of the large scale components of the energy spectrum in 
LES of hydrodynamic turbulence, this cannot be done in this 
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FIG. 1: x(0 averaged in time for jz in forced MHD turbulence. The 

^ MHD simulation, diamonds to the 
LAMHD run. The 
dashed line indicates a slope of 0.50. The arrows indicate the inertial 
range. Note that the slopes are of import, not the offsets. 



pluses correspond to the 1024 
512^ LAMHD run, and triangles to the 256^ 



context since wide energy spectra and large amounts of spa- 
tial statistics are needed to properly compute the cancellation 
exponent (see e.g. fl^ for a study of intermittency in LES). 

Fig. n shows for the three simulations, averaged us- 
ing 1 1 snapshots of the current covering a total time span of 
20 turnover times in the turbulent steady state. A power law 
can be identified at intermediate scales, scales smaller than 
the forcing band but larger than the dissipation scale. Note 
that the two LAMHD simulations reproduce the same scaling 
as the MHD simulation. As a result, the sign singularity and 
fractal structure are both well captured in the inertial range al- 
though the alpha-model is known to give thicker structures at 
scales smaller than a due to the introduction of the smoothing 
length L2i. il7ll . The best fit for the current jz using a power law 
in the inertial range gives k = 0.50±0.17for the 1024^ MHD 
run, K = 0.55 ± 0.19 for the 512^ LAMHD simulation, and 
K = 0.55 ± 0.43 for the 256^ LAMHD simulation. Note that a 
value of K = 0.50 in the MHD simulation gives a value of the 
fractal dimension D = 1.00 ± 0.34, close to the codimension 
of 1 corresponding to current sheets in MHD turbulence. For 
the vorticity, the cancellation exponent is k = 0.73 ± 0.16 for 
the 10242 MHD run, k = 0.74 ± 0.32 for the 512^ LAMHD 
simulation, and k = 0.80 ± 0.32 for the 256^ LAMHD sim- 
ulation, giving a fractal dimension of D = 0.54 in the MHD 
simulation. The values obtained are compatible with the val- 
ues of K 0.43 ± 0.06 and £» = 1.14 ± 0.12 for the current, 
and K = 0.69 ± 0.12 and D = 0.62 ± 0.24 for the vortic- 
ity obtained in Ref. fT3l for forced direct numerical simula- 
tions of 2D MHD turbulence using a 1024^ spatial grid and 
?/ = i^ = 8xl0^''. Given the good agreement between MHD 
and LAMHD simulations, in the following we will only refer 
to the cancellation exponent for the current density. 

Fig. |2]a shows the coiTesponding results for free decay- 
ing MHD turbulence. Three simulations are shown, one 
MHD run using 2048^ grid points, a 1024^ LAMHD run with 
a « 0.0058, and a 512^ LAMHD run with a w 0.0117. The 
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FIG. 2: (a) x(0 at t = 4 in the free decaying simulations, pluses 
correspond to the 2048^ MHD simulation, diamonds to the 1024^ 
LAMHD run, and triangles to the 512^ LAMHD run (the dashed line 
indicates a slope of 0.52 and the arrows indicate the inertial range); 
(b) time history of the cancellation exponent (thick lines) for the three 
runs, and of rj ( j?), where the brackets denote spatial average. 



three simulations were started with the same initial conditions; 
initial velocity and magnetic fields with random phases be- 
tween k ~ 1 and k = 3 in Fourier space, and unit r.m.s. 
values. The kinematic viscosity and magnetic diffusivity used 
were v = rj = lO^**. The three simulations were evolved in 
time without external forces. 

The evolution of the cancellation exponent as a function of 
time in the free decaying simulations is shown in Fig. |2b. 
For these simulations, the cancellation exponent is computed 
between the lengthscales L/l w 20 and L/l « 70, where a 
power law scaling in x(0 ^e clearly identified from t = 
2.5 up to t = 10. At t ~ the cancellation exponent k is 
zero, which corresponds to the smooth initial conditions. A 
gap between t = and t = 2.5 is present where no clear 
scaling is observed. As time evolves, k grows up to 0.75 at i « 
8, as the system evolves from the initially smooth fields to a 
turbulent state with strong and localized current sheets. After 
this maximum, the exponent k decays slowly in time. The 
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FIG. 3: Time history of k, (solid line) and 77 (jf) (dotted line), for a 
free decaying LAMHD simulation with r) = u — 2 x 10~^. 



maximum of k takes place slightly later than the maximum 
of magnetic dissipation, as is also shown in Fig. |2]b. Note 
that the alpha-model also captures the time evolution of the 
cancellation exponent in free decaying turbulence, as well as 
the fractal structure of the problem as time evolves. 

As previously noted in the alpha-model slightly over- 
estimates the magnetic dissipation. Note however that in the 
three simulations the peak of magnetic dissipation takes place 
close to i « 6, just before the peak of the cancellation expo- 
nent K. From the maximum energy dissipation rate, the Kol- 
mogorov's dissipation wavenumber for the kinetic and mag- 
netic energy at t « 6 are estimated as A:^ « fc^ w 470, and 
this is again larger than the largest wavenumbers resolved in 
the two LAMHD simulations. 

The observed slow decay of the cancellation exponent 
(compared with the square current) is related to the persis- 
tence of strong current sheets in the system for long times, 
even after the peak of magnetic dissipation. The system, in- 
stead of evolving fast to a smooth solution at every point in 
space, keeps dissipating energy in a few thin localized struc- 
tures. The existence of these current sheets at late times can 
be more easily verified in simulations with smaller viscosity 
1/ and diffusivity 77. While in the peak of magnetic dissipation 
the system is permeated by a large number of small current 
sheets, at late times only a few current sheets are observed 
isolated by large regions where the fields are smooth. 

Given the good agreement between direct numerical simu- 
lations (DNS) and LAMHD as seen in the preceding figure, 
we can reliably explore with the model Reynolds numbers 
unattainable in a reasonable time with DNS. In this context, 
we show that the maximum values of k obtained in the simu- 
lations seem to be insensitive to the Reynolds numbers within 
a given method (MHD or LAMHD) once a turbulent state is 
reached. As an example, in Fig. |3lwe give the time history 
of the cancellation exponent and the square current for a free 



decaying LAMHD simulation with r/ 



2 X 10"^ up to 



t = 20. The initial conditions are the same as in the previously 
discussed simulations, and a w 0.0033. It is worth noting that 
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FIG. 4: x(0 at f = 3 (dots), and t = 20 (pluses), for the free 
decaying LAMHD simulation with r) = v — 2 x 10~^. The arrows 
indicate the inertial range. 



the time evolution of the magnetic dissipation in both decay- 
ing runs (Figs. |2]b and|3) confirm previous results at lower 
Reynolds numbers |18, 19]: namely that the peak dissipation 
{t ^ 7) is lower for higher Reynolds numbers, while for later 
times it is quite independent of the Reynolds values. 

Fig. |4]shows x{l) for early and late times in the same sim- 
ulation. At small scales, the slope of x always goes to zero, as 
can be expected since close to the dissipation lengthscale the 
fields are expected to be smooth. However, note that as time 
evolves the scaling of x with / drifts to smaller scales, and at 
< = 20 a scaling can be observed up to Z/L « 0.005. By 
virtue of the model the scaling is wider and the slope goes to 
zero faster than in the DNS due to the larger Reynolds number. 



The statistics of sign cancellation in magnetofluid turbu- 
lence are related with intermittency and anomalous scaling 
of structure functions |!1 1], inherently associated with the dy- 
namics of the small scales, and as a result harder to model 
in truncations or closures of the MHD equations. For exam- 
ple, two-point closures of turbulence behave smoothly (they 
also have no information about physical structures since they 
deal only with energy spectra). The intermittency of LES is 
an open topic, in particular because for neutral fluids the need 
to study the three dimensional case implied until recently that 
only low-resolution studies could be accomplished for which 
intermittent structures were barely resolved (see jlA for a 
recent study). From that point of view, the present study 
in two dimensions allows for higher Reynolds number stud- 
ies. In MHD turbulence, the energy cascade being to small 
scales both in two and three space dimensions, it is hoped 
that the information gained here will carry on to the three- 
dimensional case. The new result stemming from this study 
is that the LAMHD alpha-model, although it alters the small 
scales through filtering, it nevertheless preserves some statisti- 
cal information concerning the small scales. It is able to repro- 
duce the scaling observed in forced MHD turbulence, as well 
as the time evolution of the cancellation exponent in free de- 
caying simulations and as such, it represents a valuable model 
for studies of MHD flows for example at low magnetic Prandtl 
number v/rj as encountered in the liquid core of the Earth or 
in the solar convection zone (see e.g., IJI). 
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